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Abstract. This work is a continuation of a previous pub- 
lication Thieberger & Spiegel concerning the use of two 
dimensional catalogues in unraveling the large scale dis- 
tribution of galaxies. Since the data in two-dimensional 
catalogues are more plentiful than those in three dimen- 
sions, the luminosity method can be helpful in using two 
dimensional data to unravel the nature of the three di- 
mensional distribution. In this work we analyse the Cfa 
catalogue and demonstrate that the results indicate a dis- 
tribution of two dimensions in the nearer range and three 
dimensions in the farther range. The possibility of a dif- 
ferent regime for very short distances is discussed. 
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1. Introduction 

Standard cosmology is based on the assumption that the 
universe is spatially homogeneous, at least on scales suffi- 
ciently large to justify its approximation by a Friedmann- 
Lemaitre-Robertson- Walker (FLRW) model. The high 
isotropy measured in the cosmic microwave background 
radiation (CMBR) is usually taken as strong evidence in 
support of this hypothesis. 

The structures seen in galaxy catalogues - groups, clus- 
ters and superclusters, distributed along voids, filaments 
and walls - are not viewed as contradicting this Principle, 
as the common opinion is that the scales on which the uni- 
verse is assumed to be homogeneous are much larger than 
those subtended by these structures. However, the consen- 
sus on a homogeneous feature of structures, even on very 
large scales, has never been complete (see e.g. Pietronero 
,1987). 

In more recent times, three-dimensional galaxy cata- 
logues have been made available, they support the findings 
of approximate self-similarity in the large-scale distribu- 
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tion of galaxies within prescribed scale intervals (Proven- 
zale, 1991). Since the data in two-dimensional catalogues 
are more plentiful than those in three dimensions and are 
likely to remain so for some time (Cress, et al., 1996), it 
seems worthwhile to discuss further the ways in which the 
two-dimensional data may be used to unravel the nature 
of the three-dimensional distribution. The basic issue in 
analyzing a two-dimensional catalogue goes back to work 
on stellar statistics (Trumpler and Weaver, 1953). The sig- 
nificant point in the stellar studies was the effect of the 
luminosity distribution of the observed objects on the pro- 
jected distributions. 

The general outline of this issue, which is the main 
concern of our present attempt has been given previously 
(Thieberger and Spiegel). Just as in the abovementioned 
publication, we shall leave out of the account some other 
features of the galaxy case that remain to be considered 
in detail, namely source evolution and cosmological ef- 
fects (Spedalere and Schucking, 1980; Ribeiro, 1995 and 
Celerier and Thieberger, 2000). 

In this study we limit ourselves to a three dimensional 
catalogue so as to be able to evaluate the results. There- 
fore we chose the cfa catalogue, which has been widely 
studied, and therefore where we can perform all necessary 
comparisons so as to understand more precisely our ap- 
proach. 

One of our motivations for the present study comes 
from our interest in trying to distinguish the different 
scaling regimes that may exist in the galaxy distribution, 
though on different length scales. It has been suggested 
(Murante et al., 1998) that three scaling regimes may be 
discerned in the distribution of galaxies: on the small- 
est scales the results are consistent with a distribution 
of density singularities; in the intermediate range, there 
is scaling behaviour suggestive of flat structures such as 
Zeldovitch (1970) favored; on the largest scales, the data 
indicate a homogeneous distribution of galaxies with non- 
fractal behavior. This categorization is based on the re- 
sults of an analysis of three-dimensional catalogues. How- 
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ever, it has to be admitted that those data may as yet 
not be adequate to clearly decide such issues and that 
is why we wish to consider whether we may reliably use 
two-dimensional catalogues for immediate needs. Because 
of statistical problems, we find it at the moment difficult 
to study the close region behaviour, therefore we will limit 
ourselves in this study to the two other regions namely, 
the two dimensional region and the three dimensional one. 

Before describing our calculations, we devote the two 
next sections to a short description of the correlation inte- 
gral method and to a brief outline of the influence of lumi- 
nosity effects by imposing a Schechter luminosity function 
on the objects of study and attempt to see how this influ- 
ences the distribution of the projected objects, or rather 
of their apparent magnitudes. Then we devote a section to 
a model catalogue so as to help us understand the results 
of the observational data. 

2. The correlation Integral method. 

Calculating dimensions for galaxies enters into difficulty 
for the cases when the distribution is fractal (Coleman 
and Pietronero 1992). A robust characterization of the 
structures of point sets is given by the correlation integral 
(Grassberger and Procaccia, 1983), defined as: 

where G is the Heaviside function. The inner summation 
is over the whole set of — 1 galaxies with coordinates 
Xj, i ^ i, and the outer summation is over a subset of N' 
galaxies, taken as centers, with coordinates X^. By taking 
only the inner N' galaxies as centers we allow for the effect 
of the finiteness of the sample (see e.g. Sylos Labini et al., 
1998). 

We may interpret C2(r) as N{r)/N where M{r) is the 
average number of galaxies within a distance r of a typical 
galaxy in the set. As r goes to zero, C2 should go to zero 
and, for general distributions, we express this conclusion 
as C2 oc r^^ . The exponent D2 is called the correlation di- 
mension and it is necessarily < 3 for an embedding space 
of dimension three. When D2 is not an integer, the distri- 
bution is called fractal (Mandelbrot, 1982). 

In Eq.(|^) we have a double summation. In certain 
cases, for example when we have a pencil beam, this in- 
ner region has to be very small and so sometimes one 
has to use just the distances from a single point. This of 
course results in a strong deterioration of the statistics. 
For comparison, we will perform such calculations and we 
will call it the "pencil beam method" (the dimension ob- 
tained will be denoted by I?pb, comparing it to the usual 
"Grassberger Procaccia method" (with dimension Dap)- 



3. Luminosity effects. 

In considering the distribution on the celestial sphere of 
galaxies or of radio sources the role of the intrinsic lu- 
minosities of the observed objects must be allowed for. 
Such effects appear in A^(> /), the number of galaxies 
detected whose flux density exceeds a given value, / (Pee- 
bles, 1993). As is done in Thieberger and Spiegel, we as- 
sume that the luminosity function of the objects in our 
study is independent of position and of time. 

A galaxy at distance r is bright enough to be seen at 
flux density,/, if its luminosity exceeds L. Using a charac- 
teristic luminosity L*, for normalization one can express 
the mean number of galaxies per steradian brighter than 
/ by (Peebles, 1993): 

N{>f)^ const. r^-^dri}{ATir^f/L^) (2) 

D being the dimension of the galaxy distribution, whether 
fractal or regular. ip{x), is proportional to the probability 
that a galaxy at distance r is included in the catalogue. 
For a Schechter model (Schechter, 1976), 

poo 

2P{x) = / exp-*" dw (3) 

Here a = —1.07 and x = L/L^. 

Lets consider our specific problem. We assume that 
(see Murante et al., 1998) at short distances we have a di- 
mension D=2 ("Zeldovitch panckakes") and at large dis- 
tances we have homogenous distribution (i.e. D=3). In 
this picture we have ignored the very close region where 
we probably have D=l (the singularity picture). Then we 
can modify Eq.(||), and write: 

N{> f) - / yd#(y') + r' " / y'dy4>{y^) (4) 



Here y = 47rr^//L*, and Tq = ro\/4:Trf / L^. ro is the value 
at which we change the dimension from two to three. It 
is easy to see that for the limit of small f the value of the 
slope logA(> /) vs. log / tends to D = 3, and for large 
/ it tends to D — 2. We used this equation to evaluate 
numerically the dimension, for the whole range of /, by 
calculating an array of values of / vs. A''(> /), making 
a numerical integration of Eq.(^). We obtained then the 
dimension by a least squares fit. This analysis resulted in 
D = 2.9 for very small /, and going fast towards the value 
of D = 2. , for large /. These results are of course quite 
obvious, the only point we wanted to clarify is the exten- 
sion of the transition zone between 2. and 3. dimensions. 
We found that this transition zone is larger than the zone 
where f tends to zero and that quite small values of / al- 
ready give a dimension D = 3., i.e. are already behaving 
as very large /. In the next section we will perform a num- 
ber of model catalogue calculations in order to study the 
relations between the results obtained by the luminosity 
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method versus the integral correlation method. We will 
also compare that to the 'beam pencil method' as was ex- 
plained in the previous section. The reason is that we wish 
to stress that in the "luminosity method" we are forced 
to consider only the luminosities relative to the observer, 
i.e. just for one galaxy. 

4. The model catalogues 

In order to understand the results which we obtain from 
an observational catalogue, we constructed a number of 
representative model catalogues. We took A^i points in a 
box to 500 (arbitrary units), and around each point on 
a plane we put N2 points at distances up to PL units, in 
each direction on the plane. To each point we attributed a 
'luminosity' obtained according to the Schechter function 
distribution (Eq. (|2|)). In most of the catalogues we consid- 
ered the cases were we had in our box around 20000 points 
(TVi X N2 minus the points which fall outside the box lim- 
its) . This value was chosen so as to be of the same order of 
magnitude as the Cfa catalogue. In this way we can exam- 
ine the reliability of the results for the real catalogue. Just 
for the sake of learning about the improvement obtained 
for larger catalogues, we checked also a few cases of a larger 
amount of points. In the following table we compare the 
results to the ones obtained by the use of correlation inte- 
grals (Grassberger and Procaccia, 1983). In the table we 
compare our results also to the pencil beam method. The 
dimensions are denoted by Dcp and DpB appropriately. 
In this table we give also the results obtained by perform- 
ing the dimension analysis via the luminosity effect. Using 
the log{f) versus log{N{> /)) results, we calculate the ap- 
propriate dimensions. The result is denoted by Dl m the 
table. 

In addition to the main purpose of the table, which is 
the comparison beween the different methods, for a vari- 
aty of cases, we also obtain an indication of the errors for 
the different models. The cases of pure two or three di- 
mensional models (the last four lines in the table) give 
good agreement and when the number of points is in- 
creased (see last line) the errors even for the luminosity 
case start to decrease considerably. The main reason for 
the larger overall errors in the top five lines is the conse- 
quence of our model where some points 'see' more three 
dimensional points then others depending on their loca- 
tion on the planes. This problem, being an artifact, should 
be less pronounced in the true catalogues. 

We wish to remark here that there are problems con- 
cerning the validity of the results both for very small / and 
for large /. For small / the finite size of our box results in 
the limit being / dependent, which is an artifact. For large 
/ only a smaller and smaller fraction of the points reach 
the "observer" , resulting in bad statistics. So we check the 
average dimension over the range of / where we thought 
the results were reliable. 



Table 1. Dimension calculation for different model cata- 
logues. 





N2 


PL 


Dgp 


Dl 


DpB 


150 


150 


50 


2.3 ± 0.1 


2.5 ± 0.1 


2.5 ±0.2 


300 


150 


50 


2.6 ± 0.1 


2.6 ± 0.1 


2.8 ±0.2 


300 


300 


50 


2.6 ± 0.0 


2.6 ± 0.1 


2.6 ±0.2 


6 


4800 


200 


2.2 ± 0.2 


2.3 ± 0.2 


2.4 ±0.3 


3700 


6 


50 


2.95 ± .05 


2.9 ± 0.1 


3.0 ±0.1 


1 


20k 


250 


2.00 ± .00 


2.00 ± .04 


2.01 ±.03 


20k 


1 


.1 


3.00 ± .00 


2.9 ± 0.1 


2.98 ±.05 


80k 


1 


.1 


3.00 ± .00 


2.9 ± 0.1 


3.00 ±.03 


320k 


1 


.1 


3.00 ± .00 


2.94 ± .03 


2.99 ±.02 



5. Observational Catalogue 

In this section we describe the dimensions that we obtain 
based on the three dimensional data of the Cfa Catalogue 
of redshifts (Huchra et al. ,1995), which is a compilation 
from several sources. A previous analysis using the stan- 
dard method, described in section 2, was performed, for 
an older catalogue, by Provenzale et al.,(1997). 

We wish to compare the results of the correletion inte- 
gral method and the luminosity method. The correlation 
integral method gives for the near galaxies D = 1.9 and 
for the far galaxies D — 2.6. The luminosity method gives 
for the mostly near (as explained in the previous section) 
D = 2.1, and for the mostly far ones: D = 2.8. 

In these cases we do not have a way to put error bars 
on the values. So to get a feel for the possible error, we at- 
tribute to each galaxy a luminosity obtained in the same 
manner as in the model catalogue, instead of using the 
values of the luminosities in the real Cfa catalogue. We 
chose nine different sets of random numbers. In this way 
we could obtain the error on the average and obtained 
an error of about ±0.2. These are encouraging results 
demonstrating that although the errors for the luminosity 
method are somewhat larger than the ones for the cor- 
relation integral method, they are in a reasonable range 
. Of course if one has a large two dimensional catalogue 
and a relatively much smaller three dimensional catalogue 
(as is at the moment the case) then the larger number of 
galaxies in the two dimensional catalogue will make the 
errors smaller then the one obtained from analysing a red 
shift catalogue. 

6. Discussion and conclusion 

Our calculations show that we obtain the same kind of re- 
sults both for the correlation integral method and the lu- 
minosity method. We tried to estimate the possible errors 
on the results and it seems to us that they are somewhat 
larger for the luminosity method. So if we have a three 
dimensional catalogue most probably the correlation inte- 
gral give a more reliable result, but as the two dimensional 
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catalogues are much larger, most probably the luminosity 
method is more reliable. 
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